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ABSTRACT 


The evolution of the UAVs leads to the creation of new configuration of multi-copters. This type of aircrafts has many advantages 
especially in civil and military applications because it can be used in high risk situations without endangering a human life and 
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inaccessible areas. In this paper, a computer simulation has been investigated to study the aerodynamic structure around a 
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configuration of multi-copters named octo-copter X8. This type of UAVs has eight rotors. Every two rotors are coaxial. The numerical 


model used is based on the Navier-Stokes equations with a k-e turbulence model. The resolution of these equations is based on the 
finite volume discretization. We have varied the length of the arm of the octo-copter and we have extracted the local characteristics 
of every aerodynamic structure using the software “Solidworks Flow Simulation” in order to determine the best octo-copter based 
on the stability criterion. 


Keywords: computer simulations, octo-copter, aerodynamic structure, CFD. 


1. INTRODUCTION 


In the last decades, academic and industrial applications have generated a big interest on the use of UAVs because of their 
capability to work in hard conditions where the surrounding environment is dangerous like in the case of civil and’especially military 
applications [1-2]. There are a wide range of UAVs. The first type is the plane using the lift power of the wing to maintain itself in 
the air [3]. Another used type of UAVs is the VTOL aircrafts and helicopter is an important configuration of it. A,helicopter is an 
aerodynamic aircraft which is designed to take off and land vertically, hover in the air and move to any direction. Necessary lifting 
force and thrust are created by one or more rotors, driven by reciprocating or jet engines. The classical,scheme of a helicopter: one 
large main rotor on top and a tail rotor to compensate the reactive moment [4-5]. The widespread type’ is the multi-rotors or the 
multi-copters. A particular case of helicopters is a flying platform with 3-6 or 4-8 and more rotors [6]. This type, in addition to its 
important advantages, presents many problems especially in the stability of flying and the lowvefficiency of its propellers. So, many 
researchers used numerical simulations to solve these problems. For example, Shi et al. [7] used numerical simulations to predict the 
aerodynamic flow around helicopter rotors in hover. Strawn et al. [8] presents highlights from CFD modeling and development of 
aircrafts. Caradonna et al. [9] presented the experimental and analytical studies of a model helicopter rotor in hover. Ansari et al. [10] 
cited new approaches from previous literature, classified them into many»types like steady-state, quasi-steady and fully unsteady 
methods used in the aerodynamic modeling for micro-air vehicles. From the mentioned references, it appears clearly that the design 
of the aircraft, the geometrical parameters and the external conditions have a big effect on the aerodynamic behavior of the aircraft. 
In this paper, we are interested on the study of the aerodynamic structure around an octo-copter. 


2. GEOMETRICAL ARRANGEMENT 


Figure 1 presents the chosen geometrical model. It is anyocto-copter drone with coaxial propellers. Each two propellers are 
mounted on a motor to facilitate the design. We choose to,work with X5C-02 main blades propeller because it is the type of rotors 
available in our laboratory. Figure 1.a shows the geometrical arrangements of the model, defined by a diameter of the propeller 
d=135 mm, a length of the arm L=135 mm, L=210 mm and L=285 mm, a length of the support Is and a radius of the support R. 

Figure 1.6 presents its exploded view. The white propeller is a Clockwise CW propeller and the red propeller is a Counter 
Clockwise CCW propeller. Considering, that the coaxial octo-copter can fly if each two opposite propellers rotate at the same 
direction, each two adjacent*propellers should rotate at different direction and the upper of the octo-copter should have the 
opposite configuration of the bottom The propellers are also fixed considering their availability in the laboratory. So, we have 
studied the effect of the variation of.the arm's length. In fact, we have changed the prototype arm-length to obtain the most 
accurate configuration which met the goals. The sizes of the arm’s length are equal to L=135 mm, L=210 mm and L=285 mm. Figure 
2 shows the 3D presentation)in each case of study. 


(a) Geometrical arrangements of the model 
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(b) Exploded view of the geometrical model 


Figure 1 Geometrical model 


(a) L= 135 mm 


(c) L = 285 mm 


Figure 2 3D case study 


(b) L = 210 mm 
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3. NUMERICAL MODEL 

To study the aerodynamic structure around an octo-copter, a computer simulation has been done. The used software is “Solidworks 
Flow Simulation”. This software is based on the resolution of the Navier-Stokes equations with a k-e turbulence model [11-12]. The 
resolution of these equations is based on the finite volume discretization [13-14]. The spatial discretization is obtained by following 
a procedure for tetrahedral interpolation scheme. For the temporal discretization, the implicit formulation is adopted. 


3.1. Control volume and boundary conditions 

Figure 3 shows a presentation of the control volume depicted with the boundary conditions. The pressure is equal to p=101325,Pa 
and the air molecular mass is equal to 0.029 kg/mol. Each one of the three octo-copters rotates in a rotational speed equal to 
Q=1711 rpm. The boundary conditions are the presented in the form of the wall limiting the control volume. This model gives us the 
most accurate values of force and local characteristics. 


Figure 3 Control volume 


3.2. Mesh resolution 

In this section, we have created a rotating flow simulation study by selecting a new project using wizard, adjusting the domain size, 
setting initial and local initial mesh, creating a rotating region and setting up the goals. In order to carry out the numerical 
simulations, it is important to set the mesh in every case of study. Figure 4 shows the chosen mesh types for the three cases. We 
fixed an unstructured mesh having the characteristics shown in table 1. 
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(a) L=135 mm 


(b) L=210.mm 


(c) L=285 mm 
Figure 4 Different meshes used in the numerical models 
Table 1 Selected mesh 
Length Number of fluid cells Number of partial cells 
L= 135mm 10646 4368 
L= 210mm 9351 3738 
L = 285 mm 10002 4573 


3.3. Mathematical formulation 


The model is/based on the Navier-Stokes equations [11] presented by the following formulation: 


OP 4 (Pu) _ 9 (1) 
Ot Ox, 


L 


A 
A(pu;)  O(Pum;) dp =— @ 47 ®)+8 

or oe ee ee I 
i i j 


Where: 
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P is the density (kgm), t is the time (s), x;and x; are the Cartesian coordinates, u; and uj are the velocity components (ms!) 


respectively on the i andj directions, p is the pressure (Pa), S; is the mass-distributed external force per unit mass (kg m*s%), T;; is 


the viscous shear stress tensor (Pa). 


In the case of Newtonian fluids, the viscous shear stress sensor is defined as: 
Ou. u, 2 
T= L( car + J = Ou, ) (3) 
ij Ox ij 
Ox, Ox, 


[I is the viscosity (Pa s) and iF is the Kronecker delta function. 


Following at assumption, the Reynolds-stress tensor has the form: 
Ou, Ou; 2 Ou,. 2 
R i = = 
T*=p( ‘+ 7-76 *)—* pkd (4) 
ij toe Ox 3 4 By 3 ij 
i k 


U 


U =f C, pk? () 
t u € 
20.5 
f = [1- exp(—0.025Ry) 21+ ___) (6) 
R, 
k2 
R = pK (7) 
Le 
Re pb (8) 
U 


Where [, is defined using two basic turbulence properties: the turbulence kinetic energy k (J.kg"') and the turbulent dissipation € 


(W kg”), ff is the turbulent viscosity.factor.and y is the distance from the wall. 


The two additional transport equations used to describe the turbulent kinetic energy and dissipation are presented as follows: 


A Us &k 
a(pk) , A(puk) “Sige *) “ y+8 9 
ao! ox @ ex Oo Ox 
i i k j 
0 o 
(oe) louse © quit) “a4 (0) 
Ot Ox Ox O° Ox . 
i i E j 
Where S; and S; are defined as: 
p OU; 
S,=T; —— pe +H :+P, (11) 
Ox. 


fi 
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E 
=C (fT rOu,+pCP)-C f pe (12) 
le li Ai, BB £2 2 k 


J 


S; 


Here Pg represents the turbulent generation due to buoyancy forces and can be written as: 


p --* 1a (13) 


Where gj is the component of gravitational acceleration in direction x, the constants os, Cy,Ci1,C:2,0; and 0; are, determined 


empirically and have these values: og= 0.9, Cy=0.09, Ci} =1.44, C:2 =1.92, 0; =1.3, Of =1. 


Cg is defined by Cg=1 when Pg>0 and 0 otherwise; 


f =1+ ase (14) 
fi 
f2 =1-exp(-R, ?) (15) 


4. NUMERICAL RESULTS 


In this section, four planes defined by x=0 mm, y=0 mm, y=112 mm and z=0 mm are considered to extract the necessary local 
characteristics. Particularly, we are interested on the distribution of’the average,velocity, the static pressure, the turbulent kinetic 
energy, the dissipation rate of the turbulent kinetic energy and the turbulent viscosity. 


4.1. Average velocity 


Figures 5, 6, 7 and 8 present the distribution of the average.velocity in the visualization planes defined by x=0 mm, y= Omm, y=112 
mm and z=0 mm. In each plane, it has been observed the effect of the variation of the arm’s length on the average velocity. 
According to these results, it has been noted that the average velocity is important near the propeller and decreases far from it in 
the three cases. Also, it has been observed thatthe high average velocity is concentrated in the upper surface of each propeller. 
When comparing the three cases, it is clear that, the more the length of the arm increases, the more the average velocity has a low 
value in the neighborhoods of the rotating volume of each rotor. The trajectories of stream lines shown in each plane present a big 
difference between the three cases. The vortex is far from the support of the drone in the case of L = 135 mm. The second and the 
third cases present more turbulent trajectories which mean less stability. 
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(a) L=135 mm 
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(b) L=210 mm 


Figure 5 Distribution of the magnitude velocity in the plane y=0 mm 


(a). L=135«mm 


(b) L=210 mm 


Figure 6 Distribution of the magnitude.velocity in the plane y=112 mm 


(a) L=135 mm 


(c) L=285 mm 


(c) L=285 mm 
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b) L=210 mm 


(c) L=285 mm 
Figure 7 Distribution of the magnitude velocity in the plane z=0 mm 


»(aL=135 mm 


(b) L=210 mm 


(c) L=285 mm 


Figure 8 Distribution of the magnitude velocity in the plane x=0 mm 
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4.2. Static pressure 

Figures 9, 10, 11 and 12 present the distribution of the static pressure in the visualization planes defined by x=0 mm, y=0 mm, 
y=112 mm and z=0 mm. In each plane, it has been observed the effect of the variation of the arm's length on the distribution of the 
static pressure. According to these results, it has been noted that the high pressure is concentrated in the rotating volume and 
decreases far from it in each rotor. Indeed, it has been observed that in the planes z=0 mm and x=0 mm, the higher pressure is 
located in the bottom of the octo-copter. However, its top is characterized by a low value of pressure. When comparing the three 
cases, it has been noted that the increment of the length of the arm increases also the zones of high pressure especially near the 
support of the drone. It has been noted that the distribution of the static pressure presented in each plane have adifference 
between the three cases. The compression zone near the support of the second and third octo-copter is higher than the one in the 
first octo-copter which means more stability. 
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Figure 9 Distribution of the static pressure in the plane y=0 mm 
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(a) L=135 mm 
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(b) L=210 mm 


Figure 10 Distribution of the static pressure in the plane y=112 mm 


(a) L=135 mm 


(b) L=210 mm 


(c) L=285 mm 


Figure 11 Distribution of the static pressure in the plane z=0 mm 


(c) L=285 mm 
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(a) L=135 mm 
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(b) L=210 mm 
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(c) L=285 mm 


Figure 12 Distribution of the static pressure in the plane x=0 mm 


4.3. Turbulent kinetic energy 

Figures 13, 14, 15 and 16 present the‘distribution of the turbulent kinetic energy in the visualization planes defined by x=0 mm, y=0 
mm, y=112 mm and’z=0 mm. In each plane, it has been observed the effect of the variation of the arm's length on the distribution 
of the turbulent.kinetic energy. According to these results, it has been noted that in the plane y=112 mm, the turbulent kinetic 
energy is concentrated in the two red propellers contrarily to the two others that have a little variation of this parameter. This fact is 
inversed in the second plane. From the comparison of the three cases, it has been noted that the variation of the turbulent kinetic 
energy depends with the variation of the length of the arm. Indeed, it has been observed a high value of this parameter in the 
second and the third cases of study. Also, the trajectories of turbulence shown in each plane present a big difference between the 
three cases. The trajectories of the turbulence present more vortex far from the support of the drone in the case of L=135 mm which 


presents more stability. 
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(b) L=210 mm 


Figure 13 Distribution of the turbulent kinetic energy in the plane y=0 mm 
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(b) L=210 mm 


Figure 14 Distribution of the turbulent kinetic energy in the plane y=112 mm 
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(¢). L=285 mm 


(a) L=135 mm 


(c) L=285 mm 
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(a) L=135 mm 


(b) L=210°mm 


(c) L=285 mm 


Figure 15 Distribution of the turbulent*kinetic energy in the plane z=0 mm 


(a) L=135 mm 
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(b) L=210 mm 
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(c) L=285 mm 


Figure 16 Distribution of the turbulent kinetic energy in the plane x=0 mm 


4.4. Dissipation rate of the turbulent kinetic energy 

Figures 17, 18, 19 and 20 present the distribution of the dissipationsrate of the turbulent kinetic energy in the visualization planes 
defined by x=0 mm, y=0 mm, y=112 mm and z=0 mm. In each plane, it has been observed the effect of the variation of the arm's 
length on the distribution of the dissipation rate of the turbulent,kinetic energy. According to these results, it has been noted that 
the variation of the dissipation rate of the turbulent kinetic energy. is located essentially at the level of the propellers. The plane z = 0 
mm shows that it is concentrated in the inferior half of the octo-copter. Indeed, it has been observed that the regions where the 
turbulent kinetic energy is high have also a high values,of theydissipation rate of the turbulent kinetic energy. When comparing the 
three cases, it has been noted that the turbulent dissipation rate depends on the length of the arm. The more the length increases, 
the more the dissipation rate of the turbulent kinetic energy increases too. When comparing the different trajectories of turbulence 
shown in the three cases, it is clear that the first case defined by L=135 mm is more stable. 
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(b) L=210 mm 


(c) L=285 mm 


Figure 17 Distribution of the dissipation rate of the turbulent kinetic energy in the plane y=0 mm 
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(b) L=210 mm 
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Figure 18 Distribution of the dissipation rate of the turbulent kinetic energy in the plane y=112 mm 
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(c) L=285 mm 
Figure 19 Distribution of the dissipation rate of the turbulent kinetic energy in the plane z=0 mm 
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b) L=210 mm 
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(c) L=285 mm 


Figure 20 Distribution of the dissipation rate of the turbulent kinetic energy in the plane x=0 mm 


4.5 Turbulent viscosity 

Figures 21, 22, 23 and 24 present the distribution of the turbulent viscosity in the visualization planes defined by x=0 mm, y=0 mm, 
y=112 mm and z=0 mm. In each plane, it has been observed the effect of the variation of the arm's length on the distribution of the 
turbulent viscosity. According to these results, it has been noted that the high value of the turbulent viscosity is concentrated in the 
rotating volume of the propeller and in his neighborhood. The plane x = 0 mm shows that the high value of the turbulent viscosity is 
situated in the bottom of the octo-copter. The more the turbulent viscosity is approached to its top, the more it decreases. When 
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making a comparison between the three cases, it has been observed that the second and the third cases do not include any 


variation of the turbulent viscosity. Which means that the increment of the length of the arm decreases the value of the turbulent 
viscosity, even eliminates it. Also, the trajectories of stream lines shown in each plane present a big difference between the three 
cases. 
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Figure 21 Distribution of the turbulent viscosity in.the plane y=0 mm 
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(b) L= 210 mm (c) L = 285 mm 


Figure 22 Distribution of the turbulent viscosity in the plane y=112 mm 
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(c) L=285 mm 


Figure 23 Distribution of the turbulent viscosity in the plane z=0 mm 
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Figure 24 Distribution of the turbulent viscosity in the plane x=0 mm 


4.6 Thrust force 

Table 2 presents a comparison between;the value of thrust force in the three cases of study and the experimental one. It has been 
noted that the more the value of the arm-length increases, the more the value of the average thrust force decreases. In our case, 
there is not a big difference between the values of thrust force in each case. When seeing the results, it has been noted that the 
experimental value of the force is near to the value of the first case of the study corresponding to the length L=135 mm. 


Table 2 Thrust-force values 


Case 1 0.081 N 
Case 2 0.067 N 
Case 3 0.056 N 
Experimental value 0.078 N 


5. CHOICE OF THE PROTOTYPE 

Figure 25 shows a 3D presentation of the final presentation of the octo-copter with the chosen length of the arm. This choice was 
made via making a compromise between the local results especially the turbulence characteristics and the value of the produced 
thrust. After this study, it has been noted that the octo-copter having an arm's length equal to L = 135 mm is the most suitable. In 
fact, according to the turbulent trajectories located in every plane, it has been noted that in the second case, the vortex is far from 
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the support of the octo-copter. As a consequence, the stability of the octo-copter will increase compared to the two other cases. It is 


true that the first case has the highest value of thrust, but it presents many zones of instability. 


Figure 25 Presentation of the final prototype 


9. CONCLUSION 

The present paper shows simulation results applied for the prediction of the thrust. As the. design of the frame is important to 
maximize the thrust and decrease the dissipation rate of the turbulent kinetic energy, the choice of the geometric model is required. 
At the beginning, we presented the proposed geometrical model of octo-copter.with the three lengths of arms. Then, we have 
extracted the characteristics of the local results which are the average velocity, the static pressure, the turbulent kinetic energy, the 
dissipation rate of the turbulent kinetic energy and the turbulent viscosity in.each case of study. According to the local results, we 
confirm that the adequate length can be determined by a compromising between the value of thrust and the stability, which means 
that the adequate frame size is the first one with L=135 mm. Indeedjit is the only case presenting vortex far from the support of the 
octo-copter which means that it is the only case that can fly with stability. 


NOMENCLATURE 

d: Diameter of the propeller (mm) 

L: Length of the arm of the octo-copter (mm) 

Is: Length of the support (mm) 

R: Radius of the support (mm) 

CW: Clockwise propeller 

CCW: Counter Clockwise propeller 

Q: Rotational speed (rpm) 

P : Density (kgm’), 

t: Time (s) 

xj and x: Cartesian coordinates 

ujand u;: Velocity. components (ms") respectively on the ( and / directions 
p : Pressure (Pa) 

Si: Mass-distributed external force per unit mass (kg ms”) 


ly : Viscous»shear stress tensor (Pa) 
L! : Viscosity (Pa s) d, is the 
d, : Kronecker delta function 


H, : Defined using k and € 


k: Turbulence kinetic energy (UJ.kg"') 
€ : Turbulent dissipation (W kg"') 


f,: Turbulent viscosity factor 
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y: Distance from the wall 


Pg: Turbulent generation due to buoyancy forces 


gi: Component of gravitational acceleration in direction x; 
Op= 0.9 

C,,=0.09 

C,, =1.44 


C5 =1,92 
0, =1.3 


O,=1 
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